Coherent transport on Apollonian networks and continuous-time quantum walks 
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We study the coherent exciton transport on Apollonian networks generated by simple iterative 
rules. The coherent exciton dynamics is modeled by continuous-time quantum walks and we cal- 
culate the transition probabilities between two nodes of the networks. We find that the transport 
depends on the initial nodes of the excitation. For networks less than the second generation the 
coherent transport shows perfect revivals when the initial excitation starts at the central node. For 
networks of higher generation, the transport only shows partial revivals. Moreover, we find that the 
excitation is most likely to be found at the initial nodes while the coherent transport to other nodes 
has a very low probability. In the long time limit, the transition probabilities show characteristic 
patterns with identical values of limiting probabilities. Finally, the dynamics of quantum transport 
are compared with the classical transport modeled by continuous-time random walks. 

PACS numbers: 05.60.Gg, 05.60.Cd, 71.35.-y, 89. 75. He, 89.75.-k 



The problem of coherent exciton transport modeled 
by quantum walks is widely studied and relevant to 
many distinct fields, such as polymer physics, solid 
statephysics, biological physics and quantum computa- 
tion P, 0, [1] . Such studies have been done in the frame- 
work of continuous-time quantum walks (CTQWs) and 
on various discrete systems 0, [1]. It has been shown 
that the dynamics of coherent transport are strongly in- 
fluenced by the structure of the underlying discrete sys- 
tems d, H • Most of previous studies focus CTQWs 
on graphs with simple structures [l^, [lH, coherent 
dynamics on general graphs have not received much at- 
tention in the scientific community. To this end, it is 
natural to consider quantum transport on graphs with 
general structure embedded in nature. 

An important and universal feature of networked sys- 
tems (or graphs) in nature is that they have the small- 
world and scale- free property [13, [S] ■ The Apollonian 
networks (ANs) [14,, i]J| are a very useful toy model that 
captures all these features simultaneously, thus provide 
a good facility to study the dynamical processes taking 
place on networked systerris, including percolation, elec- 
trical conduction, etc 14, lO. [l7j. 
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In this paper, we consider coherent exciton transport 
on 2D Apollonian networks (ANs). The network can be 
generated as follows IJ]: At the initial generation 5 = 0, 
the network is composed of three fully connected nodes 
marked as 1, 2, and 3. At the subsequent generation, a 
new node is added inside each (newly established) trian- 
gle and linked to the three vertices of the triangle. Using 
this simple rule, we can obtain a deterministic 2D ANs 
of size N = 3 + (3*^ — l)/2 (G is the number of gener- 
ation) [lj|. Many topological properties of this network 
model have been well-studied in the literature 
Fig. [1] shows the structure of an AN in four generations 
(G = 4). 

The coherent exciton transport on a connected net- 



FIG. 1: Apollonian network generated by simple iterative 
rules in four generations {G — 4). The nodes at each genera- 
tion q are marked as consecutive numbers. 



work is modeled by the continuous-time quantum walks 
(CTQWs), which is obtained by replacing the Hamilto- 
nian of the system by the classical transfer matrix, i.e., 
H = -T 1%, 20]. The transfer matrix T relates to the 
Laplace matrix by T = —jA, where for simplicity we as- 
sume the transmission rates 7 of all bonds to be equal 
and set 7 = 1 in the following ^] . The Laplace ma- 
trix A has nondiagonal elements Aij equal to —1 if nodes 
i and j are connected and otherwise. The diagonal el- 
ements An equal to degree of node i, i.e., A^ = ki. The 
states \j) endowed with the node j of the network form 
a complete, ortho-normalised basis set, which span the 
whole accessible Hilbert space. The time evolution of a 
state \j) starting at time to is given by \ j, t) ~ U{t, to)\j), 
where U{t,to) = exp[—iH{t — to)] is the quantmn me- 
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FIG. 2: (Color online) Time evolution of transition probabili- 
ties Trk,4{t) for different values of k (marked as different types 
of curves) on ANs of G = 3 (a) and G = 4 (b). The excitation 
starts at central node 4. The insets are enlarged linear-scale 
plots of return probability n4,4{t). 



FIG. 3: (Color online) Time evolution of the classical prob- 
abilities pkA{t) for different values of k (marked as different 
types of curves) on ANs of G = 3 (a) and G = 4 (b). The exci- 
tation starts at central node 4. The classical Pk,ii^) approach 
the equip-partitioned probability 1/A'' at long time scale. 



chanical time evolution operator. The transition ampli- 
tude ak,j{t) from state \ j) at time to state \k) at time 
t reads akjit) — {k\U{t,0)\j) and obeys Schrodingers 
equation [2l|. Then the classical and quantum transi- 
tion probabilities to go from the state \j) at time to 
the state |fc) at time t are given hy pk,j{t) = (k\e~*''^\j) 
and TTk,j{t) = \ak,jit)\^ - \{k\e-'*"\j)\^ 1^, respec- 
tively. Using En and |g„) to represent the nth eigen- 
value and eigenvector of H, the classical and quantum 
transition probabilities between two nodes can be writ- 
ten as [ll,|2C 



'\(ln){qn\j), 



(1) 



^,,{t) = |a,,,(i)p = |5]e-'*^"(fc|g„)(<z„|j)r (2) 



Generally, to get Pk,j{t) and ■n]^ j{t), all the eigenvalues 
En and eigenvectors are required. In the following we 
will consider Pk,j(t) and T^kjit) obtained from diagonal- 
izing the Hamiltonian H by using the standard software 
package Mathematica 5.0. 

We start our analysis by considering transport dynam- 
ics on ANs of G 3 (iV = 16) and G = 4 (iV 43) when 
the excitation starts at the central node 4. The nodes 
are numbered according to Fig. [T] and network of G = 3 
(nodes labeled as 1 ~ 16) is a subgraph of G = 4 (nodes 
labeled as 1 ~ 43). Due to rotational symmetry, the 
transition probabilities from node 4 to certain groups of 



nodes are equal. Thus, we choose several different transi- 
tion probabilities, namely, 7r4_4(i), TTi,i(t), ir^Ait), 7rg^4(t) 
and 7ri4^4(t) for further study. 

Fig-IDshows these quantum transition probabilities for 
ANs of G = 3 and G = 4. We find that there is a 
high probability to find the excitation at the initial node 
(714,4 (i) marked as solid curves in the Fig. [2]). To see 
the behavior of 7r4.4(i) clearly, we display 7r4,4(i) in an 
enlarged linear scale (See inserted plots in Fig. [2|). For 
AN of G = 3, 774,4 (i) shows regular oscillations, as gen- 
eration increases, 774,4 (i) becomes irregular and its aver- 
age value increases (Compare the inserted plots in Fig. [2] 
(a) and (b)). Transition probabilities between the initial 
node and other nodes are considerably low compared to 
the return probability 774 4 (i) (See the dashed curves in 
(a) and (b)). The corresponding transition probabilities 
7rfc,4(i) (fc 7^ 4) of a G = 3 AN is higher than those of a 
G = 4 AN. This may be attributed to the fact that the 
return probability 774,4 (i) on a G = 4 AN is larger than 
that on a G = 3 AN. 

Fig. [3] shows the classical transition probabilities 
Pk,4:{t) for different values of fc. It is found that the 
classical transition probabilities approach the equiparti- 
tion 1/7V very quickly and pi4_4(i) reaches 1/A^ much 
slower than other transition probabilities. This can be 
explained by the shortest path length from the initial ex- 
citation node 4. The number of bonds between node 4 
and 14 is larger than the distance between other pairs 
of nodes. In addition, because the shortest path lengths 
between 4 and 1, 4 and 5, 4 and 8 are equal, the classical 
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FIG. 4: (Color online) (a) Transition probabilities -Kkjit) for 
AN of G = 1. (b) Transition probabilities 71^,4 (t) for AN of 
G = 2. Both results are numerically obtained by diagonalizing 
the Hamiltonian H and consistent with the analytical results 
in Eqs. ^ and 



FIG. 5: (Color online) Time evolution of transition probabil- 
ities TTk.iit) for different values of k (marked as different types 
of curves) on ANs of G = 3 (a) and G = 4 (b). The excitation 
starts at noncentral node 1. The return probability 7ri,i(t) is 
nearly periodic for both the networks. 



Pi, 4(^^)1 ^5,4(0 Pa,4:i^) comparable (Compare the 
curves in Fig. [3]). Noting that the long time averaged 
'"'4,4 (^) is much higher than equip-partitioned probability 
and other (long time averaged) quantum transition 
probabilities is less than 1/iV, we conclude that the clas- 
sical transport to other nodes (non-initial node) of the 
network is more efficient than quantum transport. 

Interestingly, for G = 1 (iV = 4) and G = 2 (A^ = 7) 
ANs, the quantum transition probabilities are fully peri- 
odic when the coherent excitation starts from the central 
node 4. In this case we obtain, based on the analyti- 
cally determined eigenvalues and eigenvectors (2^ . that 
for G = 1, i:kj{t) have the following, periodic form: 



(5- 
(1 



3cos4i)/8, 
cos4t)/8. 



And for G = 2 we have. 



7I'fe,4(i) 



(37- 
(2- 



h 12cos7t)/49, 
2 cos 7t)/49, 



fc = 4, 



(3) 



(4) 



Fig. m shows the behavior of 'i^k,j{'t) obtained by numer- 
ically diagonalizing the Hamiltonian H for G = 1 and 
G = 2 ANs. This agrees the analytical results in Eqs. ([3]) 
and (|4]). We find that there is a perfect revival of the 
initial state for each t = 2mT/N [n e Integers), where N 
is the number of nodes of the considered network. This 
revival of the initial probability distribution resembles 
the results obtained for continuous and discrete quan- 
tum carpets 0, [l^ , in which the revival is only perfect 
for small size of cycles 0, H^l • The case for ANs is anal- 
ogous: The revivals are perfect for small ANs of G < 2, 
when the network size becomes larger (G 3), there are 
only partial revivals of the initial state (Compare Fig. [2] 
and Fig. ED. 



Now we turn to the case when the initial excitation 
starts at other positions of the network. Fig. [5] shows 
the transition probabilities when the initial excitation 
is placed at node 1. For both the G = 3 and G = 4 
ANs, the return probabilities 7ri.i(t) display regular os- 
cillations. The return probability tti i(t) is much larger 
than other transition probabilities tt^ i(t) (fc 7^ 1) at most 
time intervals. It is interesting to note that except for 
the high return probability 7ri.i(t), there is also consid- 
erable transport to nodes 2 and 3 (Note ■n2.i{t) = T^z.iit) 
because of axis-symmetry). Nevertheless, transport to 
other nodes (such as 4, 5 etc) is particularly low. This 
suggests that the excitation is preferably located on the 
nodes of the same generation of the initial node. 

If the initial excitation starts from other noncentral 
nodes, the results are similar but some details change. 
The oscillation amplitude and period are different and 
there is also a relative high probability to find the exci- 
tation at the initial node. 

In order to discuss what happens at long times, we 
consider the long time averages of the transition prob- 
abilities Pj.k{t) and 7rj_fc(i). On finite ANs, the transi- 
tion probability converges to a certain value, this value 
is determined by the long time average. Classically, the 
long time averaged transition probabilities equal to the 
equal-partitioned probability However, the quan- 

tum transport does not lead to equipartition. The long 
time average of nj^kit) is defined as 

Xfcj ^ limT^oc ^ Jq '^k,j{t)dt 

^j:n,lS{En-Ei){k\q,,){q^\j){j\q,){qi\k), ^ > 

where S{En - Ei) = 1 for En = Ei and S{En - Ei) = 
else. Some eigenvalues of H may be degenerate, so the 
sum in the equation contain terms belonging to differ- 



4 




node k 

1 ' 1 ' 1 ' 1 ' r 




node 



FIG. 6; (Color online) Long time limiting probabilities Xkj 
for different node j of initial excitation on the G = 3 (a) and 
G — 4: (b) ANs. The squares, dots, triangles and rhombus 
denote initial excitation at node 4, 1, 5 and 8 respectively. 



ent eigenstates. Here, we consider the limiting transition 
probabilities Xk.j according to this equation. 

Fig. [6] (a) and (b) show the hmiting probabihty dis- 
tributions for G = 3 and G — 4 ANs. In the figure, we 
find that Xjj is larger than other transition probabilities 
Xk,j {k 7^ j). This indicates the excitation is most likely 
to be found at the initial node, which is accord with the 
observation in Figs. [2] and [5l 

An interesting feature related to the limiting proba- 
bilities is that different nodes, k and I, may have the 
same transition probabilities, i.e., Xk,j = Xi,j- Con- 
cretely, for an excitation starting from the central node 
j = 4, transport to nodes of certain cluster has identical 
limiting transition probabilities (See the black squares in 
Fig. [6]). For instance, Xi,4: X2a and xzA are equal to each 
other; X5,4 = X&A = X7,4; Xfc,4 are equal for 8 < fc ^ 13; 
... . The nodes of clusters having the same transition 
probabilities in such case are easy to be identified due 
to rotation-symmetry of the central node. Furthermore, 
nodes of the same generation or having the same connec- 
tivity may have different limiting probabilities (compare 
the values of Xfc,4 for cluster 8 ^ /c ^ 13 and cluster 
14 < A: ^ 16 in the Fig. El). 

For an excitation starting from noncentral node, the 
situation is quite different. When the excitation starts at 
node 1 (See the dots in Fig [6]), X2,i equals to X3,i- Such 
identical values of transition probabilities are also easy 
to be distinguished and can be understood as a result of 
the axis-symmetry. The case for excitation starting at 
node 5 is analogous (See the triangles in the plots). Par- 
ticularly, if the excitation starts at node 8, nodes 10 and 



12 have the same limiting probability, i.e., Xio.s = Xi2,8 
(rhombus indicated in Fig. [6]). Such kind of identical 
probability is not straightforward to be realized but also 
can be ascribed to the rotation symmetry of the struc- 
ture of ANs. If triangle A145 is rotated 7r/3 and 27r/3, 
the initial node 8 changes to the positions 12 and 10 re- 
spectively. Except for the equal value of Xio,8 and Xi2,8 
on both the G — ?> and G = 4 ANs, we find that there 
are more identical probabilities on the G = 4 AN. For 
instance, we find the following equal transition proba- 
bilities: X31,8 = X33,8, Xl9,17 = X21,17, X26,23 = X28,23, 
Xl0,29 = Xl2,29, X31,29 = X33,29j X37,35 = X39,35j etC. For 

ANs of higher generation, the previous identical limiting 
probabilities are preserved and additional identical values 
of transition probabilities are observed due to structural 
symmetry of the network. 

In summary, we have studied coherent exciton trans- 
port modeled by continuous-time quantum walks on ANs. 
The quantum transport exhibits a very distinct behav- 
ior compared to the classical random walks. For networks 
up to the second generation the coherent transport shows 
perfect recurrences when the initial excitation starts at 
the central node • For networks of higher generation, 
the recurrence ceases to be perfect, which resembles re- 
sults for discrete quantum carpets The excitation 
depends on the initial nodes and is most likely to be 
found at the original nodes while the coherent transport 
to other nodes is particularly low. In the long time limit, 
the transition probabilities show identical values between 
different nodes, which reflects the symmetry of the net- 
work structure. 

We would like to point out that although CTQWs on 
ANs show oscillation and revivals like the results of the 
ID case, there are some difference in the quantum dy- 
namics between the two structures. For ANs, we find that 
the return probabilities at the central nodes are nearly 
periodic, in contrast to the ID case where the (maxi- 
mums of) return probability shows a power law decay as 
7r(<) - iH,!!!. In Rcf. [25], the authors find that for 
a ID chain, quantum revivals do not repeat indefin itely 
but become less and less accurate as time progresses [25| . 
They also find that the quantum walks displays Ander- 
son localizations or decoherence in the presence static or 
dynamic disorder [2^. For ANs, there are also consider- 
able localizations on the initial nodes (See Fig. [6]). Such 
localizations may relate to the network structures and 
requires a further study. 
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